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Abstract 

In this work, the dynamic behavior of the interfaces in both the standard and random 
driven lattice gas models (DLG and RDLG respectively) is investigated via numerical 
Monte Carlo simulations in two dimensions. These models consider a lattice gas of density 
p ~ 1/2 with nearest-neighbor attractive interactions between particles under the influence 
of an external driven fleld applied along one flxed direction in the case of the DLG model, 
and a randomly varying direction in the case of the RDLG model. The systems are also 
in contact with a reservoir at temperature T. Those systems undergo a second-order 
nonequihbrium phase transition between an ordered state characterized by high-density 
strips crossing the sample along the driving fleld, and a quasi-lattice gas disordered state. 
For T < Tc, the average interface width of the strips (W) was measured as a function of 
the lattice size and the anisotropic shape factor. It was found that the saturation value 
Wgg_f only depends on the lattice size parallel to the external fleld axis Ly and exhibits two 
distinct regimes: W^^^^ oc In Ly for low temperatures, that crosses over to W^^^^ oc Ly°" 
near the critical zone, aj = 1/2 being the roughness exponent of the interface. By using 
the relationship aj — 1/(1 + A/), the anisotropic exponent for the interface of the DLG 
model was estimated, giving A/ ~ 1, in agreement with the computed value for anisotropic 
bulk exponent Ab with a recently proposed theoretical approach. At the crossover region 
between both regimes, we observed indications of bulk criticality. The time evolution 
of W at Tc was also monitored and shows two growing stages: flrst one observes that 
W (xln t for several decades, and in the following times one has W (x t^' , where f3i is the 
dynamic exponent of the interface width. By using this value we estimated the dynamic 
critical exponent of the correlation length in the perpendicular direction to the external 
fleld, giving « 4, which is consistent with the dynamic exponent of the bulk critical 
transition zf in both theoretical approaches developed for the standard model. A similar 
scenario was also observed in the RDLG model, suggesting that both models may belong 
to the same universality class. 

PACS numbers:05.70.Ln; 68.35. Ct; 05.70.Np; 05.10.-a 

1 Introduction 

The study of far-from-equilibrium systems has attracted increasing interest in the last years 
[l] -[5, mainly due to the lack of theoretical frameworks, unlike their equilibrium counter- 
parts. However, useful approaches to the understanding of these systems have been developed. 
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like, for example, the analysis of simple models by means of different techniques, such as 
field-theoretical methods, numerical solving of mean-field equations, Monte Carlo simulations, 
dynamic scaling, etc. By using these procedures it is expected that one will establish and 
build more general methods to deal with these complex systems. Among the simple models 
studied, we shall focus our attention on the driven lattice gas model (DLG) introduced by 
Katz, Lebowitz and Spohn [1]. The study of the DLG model has attracted growing atten- 
tion due to its interesting far-from-equilibrium behavior. This interacting lattice gas, driven 
into nonequilibrium steady states (NESS) by an external field, exhibits remarkable properties 
such as its non-Hamiltonian nature, the violation of the fluctuation-dissipation theorem, the 
occurrence of anisotropic critical behavior [T], the existence of a unique relevant length in the 
anisotropic pattern formation at low temperatures and the consequent self- similarity in the 
system at different evolution times p], etc. 

The bidimensional DLG model is defined on the square lattice assuming a rectangular geom- 
etry (Lx X Ly) and using periodic boundary conditions. Lattice configurations are defined 
by means of occupation numbers Uij = 1 or for each site of coordinates (i, j) depending 
on whether the site is occupied or empty, respectively. By assuming nearest-neighbor (NN) 
attractive interactions between particles of the gas (J > 0) and in absence of any driven field, 
the Hamiltonian (H) corresponds to the Ising lattice gas 

H = -AJ ^ mjUi^jf, (1) 

<ij,i'j'> 

where the summation runs over NN sites only. By considering a driven field E, the coupling 
to a thermal bath at temperature T is accounted by a modified Metropolis rate, given by 

min[l,exp{AH -rjE/kbT)], (2) 

where kh is the Boltzmann constant, AH is the change in H after the attempted particle 
movement, and rj = [—1,0,1] for a particle hopping (against, orthogonally, along) the field. 
Also, the temperature T of the reservoir in contact with the sample is reported in units of 
J/hT. 

It is well known that for high enough temperatures the DLG model exhibits a lattice-gas- 
like disordered state. However, at low temperatures anisotropic NESS emerge. This ordered 
phase is characterized by the presence of strips of high particle density crossing the lattice in 
the direction parallel to the external field. So, for half-density of particles and at a well-defined 
critical bulk temperature Tc, the DLG model undergoes a second-order phase transition fl]. 
Extensive Monte Carlo simulations have shown that ~ lAlT^ , where = 2.269 J/kb is 
the Onsager critical temperature of the Ising model [1]. 

The issue of the universality class of the DLG model has become the subject of a long- 
standing debate. On the one hand, the first nonlinear Langevin equation for this model was 
developed by Janssen et ai., where it is argued that the main nonlinearity is the current term 
generated by the external field [6]. By working out this equation a set of critical exponents 
was calculated and a new universality class was found, which is different from the universality 
class of the Ising model. The values of the critical exponents were early confirmed by several 
numerical simulations [7|. Later on, another Langevin equation was developed for a variant of 
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the DLG model, called random DLG model (RDLG). In this model, the driving field is still 
applied along one axis as in the DLG model, but its direction changes randomly causing the 
net particle current, which characterizes the DLG model, to vanish. By analyzing this system 
it was found that its critical behavior falls outside of both the Ising and the DLG universality 
classes [8]. On the other hand, Garrido et al. [9] have revisited the subject by using a dif- 
ferent approach and found a new Langevin equation for the DLG model, in which the main 
nonequilibrium effect is due to the underlying anisotropy. This equation was already known 
because, in the limit of infinite driving fields (the case that commonly is studied), it describes 
the behavior of the RDLG model. According to their results, Garrido et al. have concluded 
that both DLG and RDLG models belong to the same universality class in the limit of an 
infinite driving field. These results have stimulated interesting theoretical discussions [T0| [TT]. 
and recent numerical results [El |T3] agree with the universality class of the Langevin equation 
derived from this last approach. However, subsequent numerical simulations and theoretical 
work may suggest that the DLG model belongs to the class of the Langevin equation early de- 
veloped by Janssen et al. [14] and that the RDLG model has a universality class of its own [15]. 

Besides the bulk phase transition, another interesting feature of the DLG model is the in- 
terfacial behavior and the possible existence of a roughening transition. This basically consists 
in the interface change from flat (or smooth) to rough when the temperature increases [16j . 
Originally studied in crystal growth models, it has been theoretically characterized by a Hamil- 
tonian proportional to the interface surface (the proportionality constant being the stiffness or 
surface tension), plus a periodic potential that represents the discrete nature of the crystal sur- 
face height. It has been found that the transition can be mapped to the Kosterlitz and Thouless 
transition [17j at a nontrivial roughening critical temperature, Tr, where at low temperatures 
(T < Tr), the saturated statistical width W^^^ (see Section [2]) is independent of the linear 
system size (L) parallel to it. Furthermore, for T > Tji one has W^^^f oc In L. The existence 
of the roughening transition has been confirmed in several simulations of solid-on-solid (SOS) 
models, by measuring observables such as the height-height correlation function, the difference 
of concentration in the interface layers, the fluctuations in the number of particles belonging 
to the interface [18], etc., and in experimental work concerning different metallic interfaces |17| . 

In the equilibrium Ising model, the existence of a roughening transition at a nontrivial 
Tfj depends on the system dimension. So, for d = 2 one has Tr = and for d = 3 one has 
Tr ~ 0.55 Tc (Tc is the critical temperature of the bulk phase transition) ^9\. Furthermore, 
in the rough phase (T > Tr), the structure factor G{~^) (i.e., the Fourier transform of the 
height-height correlation function {h(lt)h{0)) — (h)"^, with h{lt) being the interface height 
given at the position it), behaves as G{lf) oc l/"g^^ in the long-wavelength limit, 'qt being the 
d — I dimensional wave vector of the capillary waves [20] . This leads to a divergence in W^^^ 
for the spatial dimension d < 3. In particular, for the bidimensional Ising model it has been 
found that W^^^ ~ LP, with p = 1 [21]. 

For the case of the DLG model in d = 2, the investigations have been directed to study 
the effects of the external drive on this kind of transition. Regrettably, neither the theoretical 
approaches [22l [23] nor numerical simulations arrived at definitive results. The simulations 
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performed by Leung et al. on the DLG model at low temperatures and several values of the 
external drive E have shown that the value of the exponent p is consistent with p ^ when E 
increases, indicating the suppression of roughness, in contrast to the equilibrium Ising model 
|23| . Measurements of the structure factor of the height-height correlations are consistent with 
this picture, revealing that the singular behavior of G{lf ) is less severe than in the Ising model, 
although these data are affected by finite-size effects and it is no longer possible to give a more 
quantitative conclusion |23| . Subsequently, new measurements of G{lf) performed by Leung 
et al. using larger square lattices (the lattice size was up to 600 lattice units) showed that 
G{-t) oc l/'g'^"'', r/ = 1.33, for small t and both = 2 and ^ = 50 [24] at T ^ 2mJ/h 
<C Tc- If this behavior persists in the if —> d limit, it would be consistent with a severe 
roughness supression, due to an enhacement of the surface tension on large length scales |24| . 
In the same work, Leung et al also analyzed the interface behavior of the RDLG model in 
lattices of sizes up to 800 lattice units and at the same temperature. They found a less clear 
behavior than in the DLG model. In spite of the large size used and the noisy data near 
'g*' = 0, Leung et al. have concluded that G{lf) crosses over from G{lf) oc 1/g to G{lf) oc 
with x < 1 in the limit "g*' — > 0, suggesting again the suppression of roughness p^. From 
the theoretical point of view, some progress has been made only in the RDLG model, since 
this model has no current term in the equation of motion. It has been found that W^^-f. oc 
In L, also consistent with p = [25]. However, this dependence is not consistent with the 
simulation results for the short- wavelength limit, because a logarithmic increase of the width 
would correspond to a behavior of the type G{lf) oc 1/g in the g — > limit. 

In view of the outlined controversies and the lack of some conclusive answers to existing 
open questions, the aim of this paper is to provide an alternative evaluation of the dynamic 
critical exponents that describe the second-order phase transition of the DLG model, com- 
plementary to the one performed in ref. [I^, and to investigate the possible existence of a 
roughening transition. The attention will be focused on the properties of the interface between 
the high particle density strips and the gas phase, monitoring the behavior of observables such 
as the interface width, and the bulk and interface energies together with their fluctuations. Our 
study is based on extensive numerical Monte Carlo simulations and the interpretation of the 
obtained data is achieved by using dynamic scaling theories. We also performed simulations 
of the RDLG model in order to compare its behavior with that of the DLG model. 
This paper is organized as follows: in Section [2] the theoretical background is summarized, in 
Section [3] we give a brief description of the model; in Section [4] we present and discuss our 
results, and finally in Section [5] our conclusions are stated. 

2 Dynamic Scaling Approach for the Interface Roughness of the 
Strips 

In order to characterize the dynamic evolution of the interfaces in the DLG model, we first 
measure the average interface position (see figure [1]), given by 




(3) 
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where h{j, t) is the distance measured from the center of the strip up to the interface at the 
position of the column at time i, and Ni is the total number of particles at the interface. 
Subsequently, the interface width or roughness is defined as the rms of the mean interface 
position 



Notice that the summation in both quantities runs over all the particles of the interface because 
in general we have Ni > Ly (i.e., the number of columns given by the lattice size) and in 
consequence h{j, t) is not a single- valued function of the j-coordinate due to the presence of 
overhangs and/or bubbles. In the simulations we will compute W as a function of time, at 
different temperatures and using several lattice sizes, in order to study the change of roughness 
due to the evaporation/condensation of particles at the strip surfaces. Due to the selected initial 
condition, a single strip at T = 0, the interface has zero width at t = 0. 



Figure 1: Sketch of a strip of the DLG model (white region) showing the relevant quantities 
used for the calculations: the distance from the center of the strip (indicated by the dashed 
line) up to the interface at the position of the j^^ column at time t given by h{j, t); the average 
position of the interface at time t, {h{t)) (indicated by the dot-dashed line) (eq. jS])); and the 
average width of the interface at time t, W(t), given by eq. dl]). The gray regions represent 
the gas phase, and the external field is applied along the vertical axis. 

In the DLG model, the properties and even the existence of interfaces depend on the tem- 
perature T. In fact, for T > Tc, the system is in the disordered phase and there are no interfaces 
at all. However, for T < Tc, interfaces are present due to the strip-like patterns characteristic 
of the ordered phase (see figure [2]). If the initial configuration of the system consists of a single 
strip placed at the center of the lattice with no holes and two flat interfaces (see Section [3|) and 



W{Ly,t) = y^{hm-{h{t))^. 



(4) 
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it is taken to evolve at a temperature T > 0, the particles at the interfaces may leave the strip 
diffusing into the gas-like phase. The inverse process also takes place, and some particles of the 
gas-like phase may stick again at the interfaces. Moreover, there are two diffusion processes 
at the interface: one is driven by the external field, which sets a current along the interface, 
while the other corresponds to the diffusion of particles in the direction perpendicular to the 
driving field toward the bulk of the strip. These processes changes the interface width W (eq. 
(IH)). These changes in W come from both the intrinsic width of the interface. Win which is of 
the order of the correlation width of the bulk correlation length, and the thermal fluctuations 
of the local mean position of the interface, i.e. of the capillary waves [23]. If the system is in 
the rough phase and far from the critical point, T <C Tc, the intrinsic width is negligible and 
the interface evolution is governed by the fluctuations of the capillary waves. On the other 
hand, if T < Tc, the bulk correlation length becomes relevant, the intrinsic width is no longer 
negligible, and W oc Win. 

In order to describe the dynamic scaling properties of the interfaces of the strips of the 
DLG model, we used the scaling Ansatz early proposed by Family and Vicsek [161 [26]. So, it 
is assumed that the time evolution of roughness obeys the following scaling law 

WiLy,t)r.L^rf(^-l^j, (5) 

where f{u) oc u^' for u <C 1 and f{u) —)■ constant for u ^ 1, with Zy^ = aj/(5j. The limit 

u — > 1 sets the crossover time t^ ~ Ly between both regimes, u <C 1 and n ^> 1, respectively. 
Furthermore, a/, and z| are defined as the roughness, growth, and dynamic exponents, 
respectively. We used the subscript or superscript "I" in order to distinguish these exponents 
from those critical exponents of the second-order phase transition of the DLG model. The 
dynamic exponent describes the time evolution of the interface correlation length along the 

direction parallel to the interface according to ^| oc t^^^H [16]. Thus, for a finite sample of size 
Ly and t — > oo the correlation length remains finite ~ Ly) and the interface width reaches 
a saturation value that depends on the system size as Wsat{Ly) ~ Ly^ . Since we want to relate 
the interface behavior with the critical behavior of the system, we expect that equation ([5]) 
will hold for temperatures near the critical point (T < Tc). 
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Figure 2: Snapshot configurations corresponding to the IDLG model. Figure (a) shows a single 
strip containing all particles, corresponding to the ground state configuration (T = 0) that is 
used as initial condition in all simulations reported in this work i.e. at t=0. Figures (b)-(f) 
show the subsequent time evolution of the strip when the system is at T = 2.6 < Tc. The 
lattice size is = 100, Ly = 50 (the field is applied along the horizontal axis). Particles are 
shown in black and empty sites are left in white. Empty squares are used to show the interfaces 
of the strip. From left to right, the corresponding evolution times at which the snapshots were 
taken are the following: (a) t = 0, (b) t = 10(1, (c) t = 1024, (d) t = 10000, (e) t = 100489, 
and (f) t = 10^, measured in Monte Carlo time steps (MCS). 



In order to adapt the dynamic scaling of Family and Vicsek (eq. ^) to the DLG model, 
we have to recall that the system is intrinsically anisotropic and, consequently, the exponents 
entering into equation ^ must be consistent with that anisotropy. As we will see below, the 
width of the interface saturates for T < Tc due to finite-size effects. So, in the limit u » 1, 
equation jS]) becomes 

WsatiT) ^L^' T< T„ (6) 

where ajiT) is the roughness exponent, and Ly is the lattice size along the direction of the 
external field. 

For T = Tc and in the Ly ^ oo limit, the interface correlation length in the perpendicular 
direction of the field axis always diverges with time as a power law with dynamic exponent 
z{_ [IB]. As we will see later, the dynamic evolution of at T = Tc in large systems (see next 
Section) shows two different growing regimes (see figure [T4l) . In the first regime W diverges 
logarithmically, i.e., oc Int, while at later times W diverges as a power law, W (x t^' . So, 
according to the standard dynamic scaling theory [30], we have 



~ oc t^/^i T ~ Tc. (7) 

Since we are working with T < Tc, the major contribution to W is due to the intrinsic width, 
so Win ~ ~ Cli where is the bulk correlation length in the perpendicular direction |23| . 
Comparing equations ((H) and (jT]) for the limit u <C 1, it follows that /?/ = 
In the range < T < T^ the occurrence of a crossover time tx (see eq. ^) is justified by the 
fact that the correlation length in the parallel direction to the applied field (.^|), calculated 
from the height-height correlation function {h{'af)h{0)) — (/i)^, grows, and since the system 
size is finite in the Ly-direction, the saturation effects are observed when ^| ~ Ly for t tx- 

Recalling that oc t H we have that 

i\{tx) r^Ly(x tx'^^ . (8) 

Consequently, from the relation = a// f3j that holds for the Family- Viscek dynamic Ansatz 
[26], we have 

4 = ^ = a,4. (9) 

" PI 

From the scaling theory developed for the DLG model [Tj, it is well known that the dynamic 
exponents governing the bulk critical behavior, and z^, are related through the anisotropic 
exponent as follows [1] 

= zf/{l + AB). (10) 

We will relate these last two equations by assuming that the interface exhibits the same 
anisotropic behavior. So, the anisotropic exponent for the interface can be defined in the 
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following way: 

'" = TT^- 

The above assumptions were based on previous studies of equilibrium models where a 
conservation law is present (model B), and also in the DLG model. In fact, the bulk modes are 
slow and couple to the interface modes affecting its behavior [27l[28]. In the DLG model, it 
has been found that the bulk modes are also slow and decay slower than the interface modes, 
suggesting that they take part in the interfacial behavior [22l[29]. By virtue of these findings, 
is is expectable that the anisotropy effects induced by the driving field will be also present on 
the interface, so that = A/. In this way, by measuring both z^j_ and aj it would be possible 
to determine the critical dynamic exponent z| by using equation jO]). Furthermore, we can also 
compute the anisotropic exponent A/ by using the determined value of a/. It is worthwhile 
remarking that A^ — 1 (A^ = 2) has been predicted theoretically by the alternative coarse- 
grained description of the model developed by Garrido et al [9], (by the coarse-grained field 
theoretical equation developed by Janssen et al [HI). Consequently, the measurement of the 
interface properties of the DLG model close to criticality is expected to be helpful in clarifying 
the issue of the universality class of the DLG model. 



3 The Simulation Method 

The DLG and RDLG models are simulated in square lattices of size Lx x Ly, where L^i^y) is 
the lattice length in the perpendicular (parallel) direction to the applied drive. In order to 
investigate the behavior of the interface at low temperatures, simulations were performed by 
using lattices of sizes within the ranges 62 < L^^ < 284 and 20 < Ly < 360 lattice units (LU), 

keeping the shape factor S = Ly'^^'^^^ /L^ = 0.0786 fixed. For setting this value we employed 
the critical exponents determined by Garrido et al. [9j, i.e., u± ~ 0.63 and z/y ~ 1.22. Further- 
more, we also performed a second set of simulations keeping Ly = 20 LU fixed but changing 
Lx in the range 40 < < 4000, in order to study the dependence of the interface behavior 
on Lx- Notice that 5 is not fixed at all for these lattices. 

In the case of simulations at the critical temperature, we employed a larger lattice with 
Lx = 346 and Ly = 524 LU in order to avoid the early occurrence of finite-size effects. In all 
cases the particle density is fixed at po = 1/2, so that the DLG model exhibits a second-order 
phase transition [T]. The magnitude of the driving field is kept constant at E = 50 J for all 
simulations. This value implies, for practical purposes, that we are considering the -E — > oo 
limit (see equation So, the studied models become the infinite field DLG model (IDLG) 
and the infinite random field (IRDLG) model. The Monte Carlo time step (MCS) is defined 
such as all sites {Lx x Ly) of the lattice are selected once (on the average) in order to attempt 
the displacement of a particle to an empty neighbor site. 

The temperature is measured in units of J/kj,, kj, being the Boltzmann constant. Under these 
conditions, the bulk critical temperatures of both models have been estimated in previous work 
[DliailS], namely, T^{IDLG) ^ 3.20 and Tc{IRDLG) ^ 3.16. We monitored the time evolu- 
tion olW{t) (eq. (Ill)) over many realizations of the experiment, typically 100-3000 realizations 
according to the temperature and lattice sizes. 
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The initial condition used in order to start the simulations was the ground state configuration 
(GSC), which is the configuration that the system has at T = 0. It consists of a single strip, 
free of defects and placed at the center of the lattice (see figure [2la)). Subsequently, the time 
evolution of the interface is followed, measuring its position h with respect to the center of the 
lattice (equation ([3])) and its width W (equation Q) by identifying the particles belonging to 
the interface (this process will be described below). From the time evolution of W we measure 
the saturation value Wgat by averaging over a time interval long after the crossover time given 
by equation 

The identification of the interface is performed by employing the algorithm developed in 
ref. [H] for diffusion fronts. Let us now briefly explain how the algorithm works. By taking as 
an example the snapshot of the system exhibited in figure [21 we divide the lattice in (horizon- 
tal) rows, and suppose that we want to measure the position of the interface that is between 
the {Lx/2 + l)-th and (La;)-th rows (recall that in figure [2] L^; runs along the vertical axis). 
Then, the interface identification is performed through a process of three stages. In the first 
stage, we impose that all sites in the {Lx/2 + l)-th row are occupied with particles Q, and are 
labelled with an integer number. Then, the lattice is swept sequentially from the (La;/2 + l)-th 
row to the {Lx)-ih row, and by applying the standard Hoshen-Kopelman algorithm (HK) [30] 
to nearest neighbors occupied sites the largest cluster, which corresponds to the strip running 
along the Ly direction can be identified. In the second stage the HK algorithm is repeated but 
in the opposite way, from the {Lx)-th to {Lx/2 + l)-th rows. In this case, the largest cluster 
of empty sites, linked by both nearest- and next-nearest- neighboring sites, is determined. 
At this point, the portion of the lattice we have considered, of size Lx/2 x Ly has two clusters, 
and the interface can easily be defined as those sites of the largest cluster of occupied sites 
having at least a neighbor belonging to the largest cluster of empty sites. This is done in the 
third stage. Once the overall identification process is finished, the interface height, measured 
from Lx/2, and its width are calculated. One important feature of this method is that takes 
into account the bubbles and overhangs (they can be observed in the snapshots of figure ^ 
responsible for the intrinsic width of the interface, while previous studies neglect them by 
applying a strong coarsening of the interface, smoothing it by spreading the local density on 
neigboring sites [231 [24] . 



4 Results and Discussion 

The protocol used to register the time evolution of the interface width(s) is described first. All 
simulations were performed in the following way: we started the simulations from the GSC 
configuration and then we let the system to evolve at a temperature T < Tc for different 
lattice sizes Lx, Ly that are specified to each case. Then we monitored the time evolution of 
the roughness W of each interface, stopping the simulations when either the strip no longer 
percolated along the direction parallel to the driving field, and/or more than one strip were 
detected. Then a new simulation was performed. 

^If neccessary, all empty sites on this row will become occupied 
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In order to test the algorithm used to determine the interface position, we simulated the model 
at E = 0, i.e., the Ising model with conserved (i.e., Kawasaki) dynamics. We fixed T = 1 
and used square lattices in the range 20 < = Ly = L < 100. As it was mentioned in 
the Introduction, theoretical arguments performed by Ben Abraham et al [21j and numerically 
confirmed by Leung [22] state that the saturation width behaves as W^^^ ~ L^, with p = I. 
Figure [3] shows a plot of the saturation width W^^^ versus L. Our best least squares fit of 
the obtained data gives p = 1.07(5), which is in good agreement with the already mentioned 
expectation. This result allows us to show that the algorithm is reliable for our purposes, since 
the interface can be well characterized. 



o T=1.0 

— slope= 1.07(5) 



o 
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Figure 3: Log-log plot of the saturated width W^^^ versus L for the Ising model with conserved 
dynamics at T = 1 J/ kb. 

Now we focus our attention to the DLG model. In figures [2](a-f) a typical time evolution of 
the system from the GSC initial condition is shown. As we can observe in the figures, the strip 
does not move appreciably with respect to the center of the lattice, so cluster fluctuations due 
to the transversal periodic boundary conditions are irrelevant. This prevents miscalculations 
of the interface position. Consequently, the interface width undergoes an initial increase and 
eventually reaches a saturation value at long enough times. As expected, the saturation value 
depends on the lattice size, as shown in figure [H 
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Figure 4: Log-log plots of the interface width W versus t measured at T = 2.60 < Tc for the 
IDLG model. The lattice sizes employed were (a) = 62, Ly = 20, (b) = 100, Ly = 50, 
and (c) Lx = 208, Ly = 200. The curves were obtained by averaging over 2000, 300, and 60 
different realizations in (a), (b), and (c), respectively. 

After determining the saturation value of the interface width at different temperatures and 
using lattices of different sizes, we fitted the data of W^^^ with the aid of the following Ansatz 

Wl, = AilnLy + A2Ll''', (12) 

where the coefficients A'-s are temperature dependent, and the roughness exponent is taken as 
a J = 1/3 or a/ = 1/2, according to the values of determined by Janssen et al or by Garrido 
et al, respectively (see eq. (fTTI) and the paragraph below it). The plots of the Wsat data and 
the corresponding fits using equation (fT2]) are shown in figure [5] (a) for the case aj = 1/2. 
Figured] (a) shows that the Ansatz (eq. (fT2]l ) fits the saturation values at all temperatures, 
although the fit is clearly better for lower temperatures, e.g. T < 2.96. This feature is also 
present if we take aj = 1/3 as roughness exponent. Furthermore, figure [5] (b) shows that 
the saturated width does not depend on either the perpendicular lattice size Lx or the shape 
factor 5, since the values of Wsat obtained by keeping S fixed are also included in the plot, 
and practically no difference is noticed from those corresponding to the nonfixed S values. 
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Therefore, we can conclude that our results are not affected by any particular choice of the 
shape factor. This plot also shows that the interaction between interfaces is negligible since no 
noticeable changes in W^at are observed when the distance between both interfaces of the single 
strip is increased due to the change of by keeping the density of particles po = 1/2 constant. 



13 



40 ■ 



30 ■ 



^ 20 - 





* 


T=2.05 


+ 


T=2.30 


• 


T=2.60 


♦ 


T=2.80 


► 


T=2.90 


X 


T=2.96 


□ 


T=3.10 


O 


T=3.12 


A 


T=3.13 


> 


T=3.14 


V 


T=3.15 




1.75 



1.5 



1.25 - 



• 


T=2.60 


□ 


T=3.00 


<> 


T=3.12 


> 


T=3.14 



(b) 



64 



128 



256 



512 



1024 2048 4096 



L 



Figure 5: Saturation width versus lattice sizes, (a) Linear-log plots of the saturation width 
VFg^^j versus Ly obtained for the listed temperatures, including T = 2.05 J/fej, (see the Intro- 



duction). The full lines correspond to regression fits of the measured data by using equation 
(fT2]) taking ct/ = 1/2; (b) linear-linear plots of Wsat versus and fixed Ly = 20, for the listed 
temperatures. The values of W'^^^ (figure 4 (a)) obtained with S fixed are also included. 
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The fitting coefficients of equation (fT2]) will help us determine which behavior (i.e., loga- 
rithmic or power law) becomes relevant along the investigated temperature interval. Figure [6] 
show plots of Ai and A2 as a function of T. By comparing figures [6] (a) and (b), we can observe 
that in the interval 2.05 < T < 2.96, Ai increases monotonically assuming almost the same 
value for both cases of aj. A2 is practically zero within this regime for both cases, indicating 
that the logarithmic finite-size behavior of Wgat dominates over the power law. On the other 
hand, near the bulk critical temperature (T > 3.0), Ai still grows monotonically and then 
reaches a saturation value around Ai ~ 1.75 for the case a/ = 1/2, but whereas it fluctuates 
with a decreasing trend if a/ = 1/3 is taken. However, A2 increases for both tested values of 
a, revealing that the power-law behavior becomes relevant for this temperature interval. 




Figure 6: (a) and (b) show linear-linear plots of the coefficients Ai and A2 versus the temper- 
ature, for both cases aj = 1/3 and aj = 1/2, as indicated in the legend. Data obtained from 
the fits of the data shown in figure [5] according to eq. (fT^ . 



Since both regimes are clearly distinguishable, we investigated them separatedly. In this way, 
figure [7] shows plots of W^^^ oc In Ly as obtained within the low T range. 
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Figure 7: Linear-log plot of the saturation values of the roughness W^^^ versus Ly, for very low 
temperatures, as indicated. The data are fitted by the logarithmic term in eq. (fT2]) . indicated 
in the plot by the full lines. 

The obtained results show that the behavior W'^ ~ In Ly holds for all temperatures reported, 
including T = 2.05 J/k^, the working temperature used by Leung et al in their simulations 
[24] . This result had already been predicted theoretically by Zia et al. for the RDLG model 
[25], and tested by Leung et al. [24]. However, the best knowledge of ours, there are neither 
theoretical predictions nor numerical tests for the IDLG model that we can use for the sake of 
comparison. These results are consistent with a roughening temperature for the IDLG model 
such that Tr < 2.05. 

Within the high-temperature interval, the logarithmic growth crosses over to a power-law 
behavior (see the power-law term of equation (fT2]) ) according to equation ([6|). Consequently, 
an effective roughness exponent deff can be determined for each temperature. Figure [8] shows 
log-log plots of Wg^^ versus Ly obtained for temperatures close to Tc. 
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Figure 8: Log-log plots of the saturation width Wg^^^ versus Ly corresponding to T = 3.12 (0), 
T = 3.13 (A), T = 3.14 (>), and T = 3.15 (v), respectively. The insets show log-log plots of 
Wgg^^ corrected by the logarithmic term in equation (fT2]l versus Ly (the error bars are excluded 
for the sake of clarity). The upper inset shows the corrections and power law fits if a/ = 1/3 
is taken in equation ([6]), while the lower inset shows the same for the case a/ = 1/2. In all 
plots the lines are power-law fits of the data. 

In this way, a^ff was measured for each T, and plotted as a function of the distance to the 
critical temperature. Then, we extrapolated the value of aj to Tc. Our results are displayed 
in figure [9l According to the trend observed in the plot, it is expected that a^jf — > 1/2 
when T ^ Tc, since all the determined values are over aj = 1/3 (lower dashed line in figure 
[9|), calculated by using A/ = 2 in equation (fTTI) . A linear extrapolation of the data gives 
UeffiT = Tc) = 0.4225(4), which is close to the value expected for aj if it is calculated by 
using A/ = 1. However, the data Oeff can be improved if the logarithmic term in equation 
(fT2]) is substracted from W^^^. By fitting the corrected values of W^^^ with equation ^ (see 
insets of figure E]) , we observe that all values are, again, well above a/ = 1/3 (indicated in the 
plot by the lower dashed line) calculated by using A/ = 2 (see equation ([TT]) ). The corrected 
values of Og// considering the coefficients Ai and A2 obtained by fitting the data with the 
aid of equation (fT2]l by taking aj = 1/3, are not significantly improved with respect to the 
uncorrected values. In contrast, by using the coefficients obtained for aj = 1/2, the corrected 
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values of Oeff are markedly improved and lie close to a/ = 1/2 (indicated in the plot by the 
upper dashed line). In view of these results, we can conclude that Og// approaches a/ = 1/2 
in the T — > Tc limit. This is essential for the clarification of the issue of the universality class 
of the IDLG model. By using equation (fTTl) it is possible to compute the anisotropic exponent, 
which gives A/ 1. Since we assumed in equation (fTT]l that the anisotropic relation found 
for the bulk critical behavior holds also for the interfacial behavior, we can conclude that this 
result is in agreement with the value predicted by Garrido et al on analyzing an alternative 
mesoscopic equation 
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Figure 9: Plot of the effective exponent Oe// versus \T — Td, obtained from the power-law 
fits of the data shown in the main plot of figure [8] (O)- The full line is a linear fit of the 
uncorrected data, whose extrapolated ordinate is UeffiT = Tc) = 0.4225(4) as it is shown by 
the arrow. The dashed lines indicate the predicted theoretical values for a/: the upper line 
is the value obtained by taking A/ = 1, while the lower line shows the value of aj obtained 
by taking A/ = 2. Also the values Oe// when W^^^^ is corrected by the logarithmic term in 
equation (fT2l) are included (A for aj = 1/3, □ for aj = 1/2, respectively). 

The behavior of W^^^j for T > 3 can be further investigated by studying the bulk and 



18 



interface "energies", given by equation ([T]), and their fluctuations, defined as 

{EB,lf, 



(Ell) 



(13) 



where Ebj is the bulk or interface energy (subscripts B or I, respectively). The obtained 
results are displayed in figure [TUl (a) and (b), respectively. 
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Figure 10: (a) and (b) show plots of the bulk and interface normalized NESS energies corre- 
sponding to different sizes Ly (indicated in the legend) versus the temperature T. Here 
and Ej are the energies of the initial GSC configuration. The change of convexity in the plots 
of figure 9 (a) is indicated by the arrows on the T axis. 

As the temperature increases, the bulk energy (figure HO] (a)) decreases due to the presence 
of holes in the bulk and exhibits a change in the convexity for 2.8 ^ T < 3.1, as indicated by 
the arrows in the figure. This behavior is similar to that observed in the equilibrium 2d Ising 
model. In contrast, the interface energy grows monotonically due to the interface roughening. 
On the other hand, figure [TT] (a) shows that for all system sizes parallel to the interface (Ly) 
the fluctuations of the energy as a function of T, exhibit a peak. However, the interface energy 
fluctuations still grow when both the system size and temperature increase. These phenomena, 
together with the (non-) convexity change in the (interface) bulk energy, indicate us that the 
change in the behavior of W^^^ with T is in fact a crossover to the bulk critical behavior. These 
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findings also justify the assumption that the anisotropy of the interface and bulk would be the 
same (eqs. ^ and (fTTjl l. 




T T 



Figure 11: Plots of the bulk energy (a) and interface energy (b) fluctuations. The dashed line 
in figure (a) indicates the shift of the maxima as a function of Ly. 

We also performed test simulations of the IRDLG model. In fact, figures [T2] (a) and (b) 
show plots of the saturation width versus Ly obtained for three different temperatures and 
using the same lattice sizes as in the case of the IDLG model. These figures qualitatively show 
the same behavior already observed for the IDLG model (see figured]), that is, the saturation 
width grows logarithmically with Ly at low temperatures, in agreement with the theoretical 
results of Zia et al. p5] and the numerical simulations of Leung et al. [21]. Subsequently, the 
interface width crosses over to a power-law behavior for T <Tc. 
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Figure 12: Data corresponding to the IRDLG model obtained at different temperatures, as 
indicated, (a) Linear-log plots of W^^^ versus Ly. The full lines are the fit of numerical data 
obtained by considering the logarithmic term of equation (fT2]) only, (b) Log-log plot of Wg^^ 
versus Ly. In this case, the numerical data are fitted to a power law W'^ cx Ly°'' , according to 
eq. ([6|), or to the power-law term in eq. (|12p . 

Prom the fit of the data shown in figure [12] (b), the estimated value of the roughness is 
a^ff = 0.44(2). This result is not only in agreement with the fact that aj = 1/2 since one 
has A/ = 1 in equation (fTTI) for the IDLG model, but also validates our scaling approach in 
Section [3] for the interface roughness for the IRDLG model. However, we recognize that, in 
figure [12] (b), the data have been fitted within one decade only, because either the strip used 
as initial condition no longer percolates, or eventually it splits into more than one strip. 

At this point, and before continuing with the description of our results concerning the 
dynamic critical behavior of the interface, we would like to stress that the above-summarized 
results exhibit some discrepancies with those previously obtained by using Monte Carlo simu- 
lations by Leung et al [Mj- Our results showed that for both the IDLG and IRDLG models 
the roughness grows logaritmically with the longitudinal size Ly at low temperatures T <^Tc, 
~ InLy, while near the critical point T <Tc it increases as a power-law W"^ ~ Ly' . On the 
other hand, Leung et al found that the roughness is supressed for both the DLG and RDLG 
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models at low temperatures (T <C Tc) and large lattices H. Although we cannot give a definitive 
answer to explain the origin of these differences, we can state a few possible assumptions, which 
may deserve further research in the future: 1) The coarse grain method. In order to compute 
the interface position h{j,t) (see eq. jS])), Leung et al considered a method of smoothing the 
interface by applying a coarse graining process, which consists of the distribution of the occu- 
pation number of each site (0 or 1) between its neighbors, assigning or 1/5 to each neighbor 
and to the site itself, respectively. When the algorithm is iterated, the resulting t) is a 
single-valued function of j that takes into account only the fluctuations due to the capillary 
waves and not the contribution of bubbles and overhangs. Perhaps, by using Leung's et al 
procedure, the saturated interface width becomes oversmoothed and simplified to observe a 
logarithmic increase. In fact, we still find a logarithmic increase of W^^^ as a function of Ly in 
our simulations at the same temperature and field magnitude they reported, namely T = 2.05 
J/k}j and E = 50, respectively. 2) The error in the computation of the exponent rj ~ 1.33 
of the structure factor G{q) ~ g^^"'"^ (small "f*) for the DLG model was not reported. This 
issue is important, because if 77 has a large error, one can no longer give any conclusion about 
the nature of the asymptotic behavior in the g — > limit, and consequently to distinguish a 
logarithmic or a power-law behavior of the interface width [25] . 3) The temperature range. 
Regrettably, the results of Leung et aJ for both the DLG and the RDLG models were obtained 
for a single temperature (T = 2.05 J//cfe), in contrast to the wide range of T covered by our 
simulations. By using different temperatures we are able to clearly observe the logarithmic 
dependence of the interface width, namely W^^^j = AilnLy (figure EI). Since Ai is vanishing 
when T is decreased (2.05 <T< 2.96), it would be very difficult to distinguish a logarithmic 
dependence just by measuring the interface at a single and rather low T ; 4) The finite size. 
In their work, Leung et al used lattices of maximum longitudinal size Ly = 600 lattice units. 
However, their results were not conclusive to explain the behavior in the Ly 00 limit since 
the data for the RDLG seem to exhibit a gap in the g — > limit. On the other hand, our largest 
lattice size used was Ly = 360 lattice units, which demands a lot of computational time due 
to the algorithm that identifies the interface, the long time required to reach the saturation of 
the width, and the large number of realizations that one needs to obtain good statistics. 
It is also worth mentioning that our algorithm has succesfully been tested twice under differ- 
ent circumstances: a) for the characterization of the interface of diffusion fronts, which have 
many bubbles and overhangs. In this case the numerical results are fully consistent with the 
predictions of the well accepted standard percolation theory [M]. b) For the characterization 
of the interface between magnetic domains of the Ising model (see figure [3] of the present paper 
and related discussion). 

In order to investigate the critical interfacial properties of both the IDLG and the IRDLG 
models, we also measured the time dependence of the interface width at Tc, as shown in 
figure [131 It is worth mentioning that both models exhibit the same behavior at criticality, as 
evidenced by the overlap of the numerical data shown in figure [T3J 

brief but detailed summary of these results is included in the Introduction 



22 




Figure 13: Log-log plots of the interface width W versus t for both the IDLG and IRDLG 
models, measured at their respective critical temperatures. Results obtained using lattices of 
size Lx = 346 and Ly = 524. 

Getting further inside the time evolution of the width, two well defined growth regimes can 
be distinguished. At early times the growth is logarithmic W oc ln(t), and it lasts almost three 
decades (10 - 10^ MCS) for both models, see figures [T^a) and 1151(a). On the other hand, at 
later times, a crossover to a power-law growth regime is observed, W oct/^^ , as expected from 
equation ([7]). 
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Figure 14: Growth regimes observed during the time evolution of the interface width in the 
IDLG model at criticality. (a) Logarithmic growth; (b) Power-law growth. 
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Figure 15: Growth regimes observed during the time evolution of the interface width in the 
IRDLG model at criticality. (a) Logarithmic growth; (b) Power-law growth. 
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Focusing on the power-law behavior, we can make a fit of the data to compute the growth 
exponent of the interface, which yields Pi{IDLG) = 0.245(2), so that z{{IDLG) = 4.08(4) ac- 
cording to equation jT]). On the other hand, by fitting the power-law divergence of the interface 
width for the IRDLG model, the growth exponent Pi{IRDLG) = 0.247(2) is obtained. There- 
fore, z^j_{IRDLG) = 4.05(3). These values of the exponents strongly suggest that z^j_ = for 
both models. Also, as expected, the determination of the dynamic exponent z^j_ is not helpful 
in assigning the universality class of the studied models since, within error bars, theoretical 
calculations and numerical results are consistent with z^ ~ 4, for both models and approaches. 
In any case, these results further validate our scaling approach for the interface roughness. 

5 Conclusions 

In this work we studied the dynamic and steady-state behavior of the strip-gas interfaces in 
both IDLG and IRDLG models. We measured the average width and the roughness exponent 
aj for a wide range of temperatures T < Tc, as a function of the lattice sizes and the shape 
factor. Starting with the IDLG model, the interface width reaches a saturation value W^^^ 
that depends on the lattice size parallel to the external field direction, Ly (figure [5] (a)), but it 
is independent of both the transversal size and the shape (figure [5] (b)). In the low-temperature 
regime, T < 3.0, W^^^ exhibits a logarithmic dependence on Ly, as has been theoretically found 
by Zia et al. ^5] and tested numerically in a subsequent work by Leung et al. [24] for the 
RDLG model. This behavior is consistent with the presence of a roughening transition at very 
low T, presumably at T = as in the RDLG model [25]. Then, close to the critical zone, 
3.00 = 3.20, W^^f crosses over to a power law of the form W^^^ oc L^"^^ . The value 

of Oe// is estimated to be a/ = aeff{T — > Tc) s::: 1/2 so A/ ~ 1 using equation (fTT]l . Since we 
assume that the bulk anisotropic behavior holds also for the interfacial behavior, this value is 
in accordance with that computed by Garrido et al. for the bulk anisotropic exponent A^. 
Furthermore, the field-theoretical predictions of Janssen et al for this value (A^ = 2) seems 
to be ruled out. 

In order to investigate the change of behavior in W^^^ around T = 3.00, we measured the 
bulk and interface energies, together with their respective fluctuations. The curves of the bulk 
energy versus the temperature show a change in convexity around the size-dependent critical 
temperature. The energy fluctuations are peaked at the finite-size critical temperature, and 
they shift to the — > oo critical temperature when Ly is increased. In contrast, the interface 
energy and its fluctuations grow monotonically as the temperature and the lattice size are 
increased. These facts allow us to conclude that the change in the behavior of W^^^ around 
T = 3.00 is in fact a crossover to the bulk critical behavior, and apparently it is not related to 
any interface phenomena such as a roughening transition. 

Later, the time evolution of the interface width at the critical point was measured. In this case 
W does not saturate and exhibits two growing stages. In the first stage, W grows logarithmi- 
cally for several decades, W (x In t, and in the following stage W grows obeying a power-law, 
W (xt^' , with growth exponent /?/ = l/z^_^. The obtained value of the dynamic exponent was 
z^j_ = 4.08(4), which is in good agreement with the theoretically estimated values for the bulk 
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critical behavior, that is ^ 3.996 p] and z;^ f» 4 [6]. 

On the other hand, the same studies performed for the IRDLG model gave practically the 
same results. At low temperatures Wg^^ oc InLy, confirming the predictions of Zia et al. |25] 
about the existence of a roughening transition at T = 0. Then, W^^^ crosses over to a power 
law for T <Tc, with an effective roughness exponent that is similar to the one measured in the 
IDLG model at the same temperature. At the critical point, the time evolution of W is the 
same as in the IDLG model, exhibiting the same logarithmic and power-law growth regimes. 
In this last stage the dynamic exponent of the evolution was obtained giving z^_^_ = 4.05(3), in 
agreement with the estimated value of for the IDLG model [Ull]. 

In short, we think that our study is valuable because it exploits the saturating character of 
the interface in the T <Tc interval, extracting useful information about the interface behavior 
relating it to the bulk critical behavior. Furthermore, our generalization of the dynamic scaling 
approach of Family and Vicsek to the driven lattice gas allows for an independent estimation 
of the anisotropic exponent, further estimulating the discussion about the universality class of 
the model. 
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